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Abstract 



Contributions to HFS and to the Lamb shift intervals of order 
o?'{Zaf' induced by graphs with two radiative photons inserted in the 
electron line are considered. It is demonstrated that this last gauge 
invariant set of diagrams which are capable of producing corrections of 
considered order consists of nineteen topologically different diagrams. 
Contributions both to HFS and Lamb shift induced by graphs con- 
taining one-loop electron self-energy as a subgraph and also by the 
graph containing two one-loop vertices are obtained. 
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1 Introduction 



Recent theoretical work on high order corrections to hyperfine sphtting and 
Lamb shift concentrated on calculation of nonrecoil contributions of order 
a'^{Za)^. These corrections are of immediate phenomenological interest for 
HFS measurements in the ground state of muonium and for n = 2 Lamb shift 
measurements in hydrogen. Their magnitude may run up to several kilohertz 
for both intervals to be compared with current experimental uncertainty 
equal to 0.16kHz for HFS in muonium and to 1.9kHz for the Lamb shift 
in hydrogen 0. Further experimental progress is envisaged at least in the 
case of muonium hyperfine splitting measurements [0. 

The only terms, which are capable to produce corrections with magni- 
tude about several kilohertz to muonium HFS (an order of magnitude larger 
than the experimental error) and are still uncalculated are of functional form 
a'^{Za)EF (see, e.g. For the Lamb shift in hydrogen corrections of 

order a'^{Za)^m are also as large as several kilohertz but in this case there 
are also other sources of contributions of comparable magnitude (see, e.g. 
i). 

We have shown recently that there exist six gauge invariant sets of dia- 
grams (see Fig.l), which produce corrections of order a^{Za)^ to muonium 
HFS 0. All these diagrams may be obtained with the help of different ra- 
diative insertions from the skeleton diagram, which contains two external 
photons attached to the electron line. Contributions induced by polariza- 
tion operator insertions in external photons and by simultaneous insertion of 
radiative photon in electron line and one-loop polarization operator in exter- 
nal photon had been calculated in analytic form ||^. Correction induced by 
polarization operator insertions in radiative photons was obtained in semian- 
alytic form as one-dimensional integral, where integrand is itself a complete 
elliptic integral Contribution induced by the fifth gauge invariant set 
of graphs containing light by light scattering insertions was reduced to the 
three dimensional integral over Feynman parameters which was calculated 
numerically [0. 

Contributions to the Lamb shift of order a'^{Za)^ are induced by the same 
six gauge invariant sets of diagrams which were just described for hyperfine 
structure. The only difference is that the tensor structure which is relevant 

^Several leading logarithmic corrections of higher order which nevertheless turn out to 
be numerically significant were recently obtained in 
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for the Lamb shift differs from the one which was relevant for HFS. Numerical 
calculation of contributions induced by diagrams with polarization operator 

insertions in the external photons and in the radiative photons as well as 

by the graphs containing light by light scattering insertions was performed 

recently 0, |TT] 



and m 



In the present paper we discuss all possible contributions to HFS and to 
the Lamb shift of order a'^{Za)^ which are induced by the last and most 
bulky set of diagrams with two radiative photons inserted in the electron 
line (see representative graph in Fig.lc). We present below results of calcu- 
lation of all contributions both to HFS and to the Lamb shift induced by 
diagrams containing one-loop electron self-energy as a subgraph and also by 
the diagram containing two one-loop vertices0'[|. 



2 General Strategy of Calculations 

The set of diagrams with insertions of two radiative photons in the electron 
line which is relevant for calculation of corrections of order a^(Za)^ to HFS 
and the Lamb shift contains nineteen topologically different diagrams and 
is presented in Fig. 2. The simplest way to describe these graphs is to re- 
alize that they were obtained from three graphs for the two-loop electron 
self-energy by insertion of two external photons in all possible ways. Really, 
graphs 2a — 2c are obtained from the two-loop reducible electron self-energy 
diagram, graphs 2d — 2k are the result of all possible insertions of two ex- 
ternal photons in the rainbow self-energy diagram, and diagrams 21 — 2s are 
connected with the overlapping two-loop self-energy graph. 

As was discussed at length in ^ (see also for the case of the Lamb 
shift) one has to perform a rather tedious analysis to find out which diagrams 
and in what kinematic conditions are relevant for calculation of the energy 
shifts. Happily, to this end we may use the results of just mentioned papers 
and we know that to obtain contribution to the energy splitting one has 

^Prof. T. Kinoshita announced recently in the talk presented at the International 
Workshop on Low Energy Muon Science, Santa Fe, NM, April 4 - 8, 1993 preliminary 
result of his calculation of the contribution to HFS of order a'^{Za)^ induced by graphs 
in Fig.lc (Preprint CLNS 93/1219, May 1993). 

■^Dr. K. Pachucki is also working now on calculation of contribution of order a^(Za)^ 
induced by graphs in Fig.lc to the Lamb shift. (Private communication from K. Pachucki 
to M. Eides). 
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to calculate matrix elements of the diagrams in Fig. 2 between free electron 
spinors with all external electron lines on mass-shell, project these matrix 
elements on the respective spin state and multiply the result by the square 
in the origin of the Schrodinger-Coulomb wave function. 

It should be mentioned that some diagrams in Fig. 2 contain also contri- 
butions of previous order in Za. Physical nature of these contributions is 
especially transparent in the case of HFS. As was mentioned in they cor- 
respond to anomalous magnetic moment, their true order in Za is lower than 
their apparent order and they should be subtracted from the matrix elements 
prior to calculation of the contributions to HFS. Analogous situation holds 
also in the case of the Lamb shift. The only difference is that this time not 
only the slope of the Pauli formfactor but also the slope of the Dirac form- 
factor of the electron is capable to produce lower order contribution to the 



splitting of the energy levels (see, e. g. 0, [|TT| and |T2|) and all respective 
terms have to be subtracted. Technically cases of lower order contributions 
both to HFS and to the Lamb shift are quite similar. Lower order terms are 
each time proportional to the exchanged momentum squared and to get rid 
of them one has to subtract all low frequency terms proportional to the ex- 
changed momentum squared in the asymptotes of the matrix elements when 
they exist. 

Actual calculation of matrix elements of the diagrams in Fig. 2 is impeded 
by the ultraviolet (UV) and infrared (IR) divergencies. To get rid of UV 
problems we work only with renormalized (subtracted) graphs and subgraphs 
and only such graphs are presented in Fig. 2. This is the reason why some 
obvious diagrams with e.g. self-energies on the external electron lines are 
absent in Fig. 2. IR problems are as usual more difficult to deal with than the 
UV ones. We choose the Fried- Yennie (FY) gauge for the radiative photons. 
The advantage of this choice of gauge is connected with maximally mild 
behavior of all matrix elements in the IR region of integration momenta in the 
FY gauge and with the possibility to perform on-mass-shell renormalization 
without introduction of the infrared photon mass (see,e.g. Q]). This last 
feature of the FY gauge is especially important for us as we want to push 
analytical calculations as far forward as possible and working without the 
photon mass gives us the chance to obtain much simpler analytical formulae. 
On the other hand working in the FY gauge without the photon mass makes 
it absolutely necessary to pay special attention to the infrared behavior of the 
integrand functions and to perform cancellation of spurious IR divergences 
prior to integration. 
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Each contribution of order o?{Zaf' arises from radiative insertions in the 
skeleton graph in Fig. 3 with two external photons. Contribution to hyperfine 
splitting, produced by diagrams in Fig.2 is given by the expression (see, e.g. 

i) 

A£„.,, = 8^A^J5.r^F(k). (1) 

Txn-^ Tx Jo 

where fc^ = (0, k) is the momentum of the external Coulomb photons 
measured in the electron mass units^ and 

where = m/(l + m/M) is the reduced mass of the electron-muon 
system and is the anomaly of the muon magnetic moment. The function 
F(k) is connected with the numerator structure of each particular graph and 
describes radiative corrections to the skeleton diagram. It is normalized on 
the skeleton numerator contribution. 

As was already discussed above anomalous magnetic moment of the elec- 
tron which is also produced by the radiative insertions in the electron line 
in the graphs in Fig.2 is not sufficiently mild near the lower boundary of 
integration region over exchanged momentum in eq.(0) and naively leads to 
divergence of the integral for the energy splitting. This simply means that 
anomalous magnetic moment produces contribution of lower order in Za and 
respective terms should be subtracted from the expression for the electron 
factor -F(k) prior to integration. 

Contribution to the Lamb shift induced by the diagrams in Fig.l has the 
form (see, e.g. [|^) 




Function /^(k) as in the case of the contribution to HFS is calculated for 
each particular graph and describes radiative corrections to the skeleton di- 
agram. Spin structure which is relevant for the Lamb shift differs, of course, 
from the one relevant for the HFS. The integral for the Lamb shift interval 
in eq.(^) is more singular in the infrared region than the one for the HFS 
in eq.(|I|). This singularity once again indicates that some diagrams in Fig.2 

^AU momenta below are measured in the electron mass units if reverse is not stated 
explicitly 
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contain not only contributions of order a'^{Za)^m to the Lamb shift but 
also contain abundantly contributions of lower order in parameter Za (com- 
pare the contribution of the anomalous magnetic moment of the electron to 
HFS discussed above). These contributions are connected with integration 
over external photon momenta of characteristic atomic order mZa and they 
should be subtracted to make calculation of corrections of order a'^{Za)^'m 
possible 0. 



3 Calculation of Simplest Contributions 

3.1 Diagram with Two One-Loop Self-Energy Inser- 
tions 

3.1.1 Contribution to HFS Interval 

We begin actual calculation with the simplest possible graph in Fig. 2a, con- 
taining reducible electron self-energy insertion in the skeleton graph. Explicit 
expression for the one-loop electron self-energy in the FY gauge has the form 
(see, e.g. 

E^^ip -k) = {p-k- l)^(^^^g^)M(k), (4) 

where 

^(k) = Y^+(^^logk^ (5) 

and kinematic conditions are defined by relations = (1,0), pk = 0, 
k, = {0,k). 

Substitution of this expression for the self-energy operator instead of the 
electron factor -F(k) in eq.(|l|) leads to the one-dimensional integral for the 
contribution to the HFS 

^Kfs = ^^M-7^,) rrf|k|(l-k^)M2(k), (6) 

which may be easily calculated analytically 

4 TTW^ 
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3.1.2 Contribution to the Lamb Shift 



It is not difficult to see that due to the tensor structure of the electron-hne 
factor with two one-loop electron self-energy insertions it does not lead to 
nonvanishing contribution to the Lamb shift of the order under consideration 



3.2 Diagrams with Simultaneous Insertions of One- 
Loop Electron Self-Energy and Vertex 

3.2.1 General Expression for the One-Loop Vertex with One On- 
Mass-Shell Leg in the Fried- Yennie Gauge 

Renormalized vertex operator with on-mass-shell left leg enters as a subgraph 

in the diagram in Fig.2b. We obtained general expression for those entries 
in such vertex operator in the FY gauge which produce contributions to the 
hyperfine splitting and the Lamb shift (terms which are proportional to k/^ 
are omitted because they are irrelevant for our goals) 



AEl = 0. 



(8) 



A 



FY 



(9) 



where 




(11) 



(10) 



C(k) 



l-m-s{k)-l-M{k)+hogk'- 



3 logk^ 
21 -k2 



(12) 



E{k) = -M$(k) - l5(k) - lM(k) 



(13) 
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13 2 3 1ogk^ 
+ 4 + 8^°^^ -81^ 

and 

m = iki t log = Li(i - |k|) - Li(i + |k|) (14) 







+2 



Li(l + ' ' ' ) - Li(l 



vwn, Vk^+iki 

= oil I log- 
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2|k| Vk2 + 4-|k| 

Euler dilogarithm Li is defined here as in [Q. We would hke to mention 
that function $(k) emerged for the first time in calculation of contribution 
to HPS induced by the diagrams in Fig.lc 

3.2.2 Contribution to HFS Interval 

It is not difficult now to obtain contribution induced by diagrams in Fig. 2b 
to HFS. Taking into account combinatorial coefficient 2 we obtain 

2AE'^^s = ^-^EA--,) r d\k\M{k)[-k'A{k)+B{k) + lc{k)] (15) 

7Cn'^ TT^ JO Z 

or numerically 



a\Za 



2AEljp. = -6.65997(1) ' , ' Ep. (16) 
3.2.3 Contribution to the Lamb Shift 

Contribution to the Lamb shift induced by the diagrams in Fig. 2b has the 
form 

3 ^,2/7^,\5 



'^^i - " (^) ^ (^) f^ikiMWi-^w + fl(k) (17) 

+C(k) - B(k)] 



or numerically 



2A£;^ = 2.9551(1) m — ^ / . (18) 
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3.3 Diagram with Simultaneous Insertions of Two One- 
Loop Vertices 

3.3.1 Contribution to HFS Interval 

Contribution to HFS induced by the diagram in Fig. 2c has the form 

^^^^^ = ^^^-^^^ r ^1^1 (19) 

vrn'^ TT^ Jo 
-2yl(k)^(k) - [E{\i) - B{\i)f + C(k)[£;(k) - S(k)]} 
or numerically 

A£;|fpg = 3.93208(1) °' ^^f^ Ep. (20) 

3.3.2 Contribution to the Lamb Shift 

Contribution to the Lamb shift induced by the diagrams in Fig.2c has the 
form 

+C(k) - B(k)l 

or numerically 

AE2 = -2.2231(1) mf^y^^. (22) 

4 Calculation of Contributions Induced by the 
Diagrams with Insertions of the "Left" Self- 
Energy Operator 

4.1 Diagrams Containing the Two-Loop Vertex 

4.1.1 General Expression for the "Left" Vertex with Left One- 
Loop Self-Energy Insertion 

Let us consider two-loop vertex diagram which is contained as a subgraph in 
Fig. 2d. Insertion of the one-loop electron self-energy in the one-loop vertex 
corresponds in the FY gauge to the replacement 



8 



of the left propagator in the one-loop vertex integrand. We used here one- 
dimensional integral representation for the renormalized self-energy operator 
in the FY gauge and took into account that the would be loop integration 
momentum q unlike momentum k of the external photon in eq.@ is not 
orthogonal to momentum p of the external electron. 

Formal expression for the vertex has the form 



A-K Jo 1-xJ ivr^ g2^g2 _^ 2g(p - A;) + P](g2 + 2]9g - Y^)]' 

(24) 



where 



= riP + q)l^.{P + q-k + 1)7., (25) 
Nl = 2q{P + g)7^(p + q-k + l)q. (26) 

It is convenient to put down second numerator in the form 

Nt = q'Nl, + Nt,, (27) 

where 

Nt, = 2[7^(p -k + l)q + qfn^ + q^l^l (28) 

Nt2 = '2qVl^.{v-k + l)q. 

As a result of this transformation numerator N'l^ enters all expressions 
below on the same footing as Np. Combining denominators in eq.(p^) 

(1 - t)q^ + t[u{q^ + 2g(p - k) + fc^) (29) 

+ (1 - u){q' + 2pq - -^)] = (g + Qf - A,, 
1 — X 

Q = pt — kut, 
Ai = t[k2M(l -ut)+t+ = tai, 

we obtain after shift of the integration variable q ^ q — Q representation 
for the bare vertex operator 
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^ Mtt Jo l-x Jo Jo J i 

Mi<l-Q) + Nt^{q-Q) 



{2t'- 
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Momentum integration leads to 

X f\ f\ rr^ /i A^ + 



A, 



= 3(^)2 dx^^ [' dt r du{[6taog ■ 

^ An Jo 1 — X Jo Jo 

4:t - 2t 

+—(1 - 2t + fu)pi,{p -k-l) + —[{l-t)^- ut]-f^k 
+2(1 - + ^(-Sk^iit + 2t)h, 

+^2t{l - t)[-kMi + 1]7m(p - ^ - 1) 

+4(i-t)[|^-^tK(p-^-i) 

+2(1 - ^)[^ - ^«^(k'(l - 2«) + 2)]^,k}. 
Subtraction term is equal to 



where 

A;o = taio = A;|k=o- 
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Note that we preserved term A when it is added to the ultraviolet cutoff 
in eq. (|3l|) and eq.(32) in the terms which do not vanish when the cutoff 



goes to infinity. This was necessary to preserve convergence of the subsequent 
integration over the Feynman parameter x in eq.(|3lD and eq. ([3^) in spite of 
the singular nature of the explicit integration weight x/{l — x) near x = 1. 
Really, it is easy to check that A itself is singular near x = 1 

A^i^^i ^ A^i.^i ^ ^SLlA (33) 
' I — X 

and compensates would be divergence. All potential divergences in inte- 
gration over X vanish after subtraction and subtracted vertex admits the limit 
of infinite cutoff. All extra terms containing ultraviolet cutoff A^ which were 
preserved in eq. (|3lD and eq. (|32[) to make the integration over x finite vanish 
in the subtracted expression for the vertex at least as (1/A^) log(A^). Hence, 
subtracted expression for the vertex coincides with the naive one which would 
be obtained if one simply missed the problem of finiteness of integration over 
X in eq. (pl|) and eq. (p^) ! 



Last terms in eq.(|3T|) and in eq. (^2]) may be simplified with the help of 
identity in eq.( p.3U| ), where one has to substitute n = I, r] = 1 and r = u. 
We also simplify logarithmic term in eq.(^) and eq. (|32|) which admits (as 
was just explained) omission of the ultraviolet cutoff after subtraction. We 
perform integration by parts over u, under which subtracted logarithmic term 
transforms as 

air, , ,,k^(l — 2ut) X ,\ 1 , , 

eUog ^ ^ -6t 1 - u){-^ ^ + 34 

ai ai 1- X aiQ ai 

= 6t—[l - 2ut - uH + ^]. 

ai aiQ 



Sum of the expressions in eqs. (|3^) and (0) and other terms in eq.(|3l|) 
proportional to 7^ is equal to 

Ai = - — {Q + Au- 12ut + 2uH) (35) 
6tM'^-ut) At{l - t)u{l - u) 



-i 

aiaio a I 

We also introduce several additional functions in such way as to make 
final expression look like the one for the one- loop vertex in eq.(P) 
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Bi = 2t[- + ^(1 - k\{l - u))], (36) 

CLl CLj^ 

2-^3t+fu t{l-t) 

' - n ^2 J' 

El = 2 ^ r^[k'(l - 2u) + 2]. 

ai ai 

Then renormalized vertex with self-energy insertion has the form 



dx^— C dt C du{]i^Aaa (37) 
Jo 1 — X Jo JO 



+Bimjf,{p -k-l) + CiPf,{p - A; - 1) + Sa^k}. 

Note that all four functions Ai{\<.), Bi(k), Ci(k) and Si(k) are finite at 
k = 0. 



4.1.2 Contribution to HFS Interval 

It is not difficult now to obtain contribution to HFS 

2AE^HFS = "'^y ^F(^) dx-^ C dt C du r d\k\ [^z(k) (38) 
vrn'^ vr"^ JO 1 — x Jo Jo Jo 

g/(k)-g/(0) . 

where 

_ EjW-m ^ _2t ^ _ ^ _ Miz^ (39) 

ai af 

_^ 2u{l -ut) _ ^ _ _ 4^2^(^ - t)(l - Mt)(^ ^ 



a/Oio o-i 0,10 cbio-io 

Subtraction of Si{0) in eg . (1381) corresponds to subtraction of the contri- 
bution to HFS induced by the second order anomalous magnetic moment 
which is of lower order in parameter Za. 
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After momentum integration we obtain three-dimensional integral for 
the contribution to HFS induced by both diagrams with nonsymmetric self- 
energy insertions in Fig. 2d 



2AE'^ ^ f A) f' dx^- C — C (40) 

fd fd fd 

r/l/2 /3/2 h/2_. 
\ 1 "T 3 "T 5 J) 

aio^ aio^ aio^ 



where 



= -2t(3 + 2u- 6ut + uH), (41) 
= 2m(1 - ut){2 -3t + - ut) - ut{l - t), 
fi/^^-6uH(l-t){l-ut). 
Numerically we obtain 

2AEf,p.s = -3.903(l)^^^£;f . (42) 

4.1.3 Contribution to the Lamb Shift 

Contribution to the Lamb shift has the form 

2AEi = -m C^fl^(^) r (43) 

m Tm"^ vr"^ JO I — x 

f\.. r.,u,^'(k)->c'^(o) 



ri poo 

dt du I d|k| 
Jo Jo 



k2 



where 



and 



£^(k) ^ Ai{k) - isKk) - IcKk) + ^Ei{k) (44) 



— ^ = ^ h ^ ^ 2 + St(l - At) + ut(5 + Ut - lOut)] 

k^ af aiaio 

(45) 



+ t{l - t)u{l - ut){-l + 6u- Au^){— — + 
o^afo (^1(^10 
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Subtraction in eq.(^) corresponds to subtraction of the contributions 
to the Lamb shift induced by the slope of the Dirac formfactor and by the 
anomalous magnetic moment which are of lower order in parameter Za. 

After momentum integration we obtain three-dimensional integral for the 
contribution to the Lamb shift induced by both diagrams with nonsymmetric 
self-energy insertions in Fig. 2d 

2AEf = m — f — -^(—) / dx / / , 46 

^ ^m^ Tin? 1 - x Jo y/u Jo y/T^lR ^ ' 

( 93/2 95/2 
i I" ^ FJ ' 

where 

gl/^ = -2n(l - ut){2 + 3t + but - IQut^ + 2mV) - t(l - t)u^, (47) 

= 3Mt(l - ut)[Aut{l - ut) + (1 - t)(l - 6m + Au% 
Numerically we obtain 



m^,qa^(ZQ;)^ 



2AE1 = -5.235(2)m (—Y ^ / . (48) 



4.2 Spanning Photon Diagram with the "Left" Self- 
Energy Insertion 

4.2.1 Contribution to HFS Interval 

Consider now spanning photon diagram with "left" self-energy insertion in 
Fig.2e. Unlike discussion of the vertex in the previous subsection we will 
consider here only the part of the spanning photon diagram relevant for 
hyperfine splitting. Insertion of the one-loop electron self-energy in the FY 
gauge corresponds to replacement described in eq. ( p3D in the left propagator 
in the one-loop vertex integrand. 

Formal expression for the spanning photon diagram has the form 
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q2^q2 + 2q{p - k) + fc2](g2 + 2pq - j^)iq^ + 2pq) ' 
where 

N'f = YiP + q)l^{p + q-k + 1)Mp + '1 + 1)7^, (50) 

Nl = 2q{p + q)^^{p + q-k + l)-f^{p + q + l)q. (51) 

It is convenient prior to following transformations to rewrite numerator 
of the longitudinal part in the form 

Nl = 2(g2 + 2pq)[q'^-^^{p + q-k + 1)7^ + qp-f^{p + q-k + 1)7,.] (52) 

^iq' + 2pq)[q'Nl, + Nl,]. 

Hence, denominator factor (g^ -|- 2pq) in eq.(^) cancels with the same 
explicit factor in the longitudinal numerator in eq. (B2|) and we use the identity 



1 -J. 1 

dx— = (53) 

l-x(g2 + 2pg- ^)(g2 + 2pg) ^ ^ 

^ dx 1 



1 - X (g2 + 2j9g - Y3^)2 

to get rid of the same factor in the term containing numerator Np. 

Next we combine denominators as in eq. (|29|) and after the shift of integra- 
tion variable q —>■ q — Q (see eq.(^)) obtain retaining only those numerator 
structures which are relevant for contribution to HFS 

^% = ^i-? dx- du dt (54) 

^ 8 TT JO 1 — X Jo Jo 



{q^-^lY 



J ITT'^ 



txNl,{q -Q) 3t{l-t)xNl,{q-Q) 



{q^-^if {q^-^iY 
Next we project on the spinor structure relevant for HFS splittting (see, 



e.g. |jT5| for explicit expression for the respective projector) and obtain after 



momentum integration scalar electron factor 
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F' = h' C dx-^ f'du (55) 
8 JO 1 — X JO JO ai 



+a-u 



ai 

kW(l - ut) + [-1 + t{2 - - ut)] 



{l-ut)x (1 



ai af 

Substituting electron factor in eq.(^) in eq.(|l]) one obtains contribution 
to HFS interval 

2AE^HFS = ^Q—,{-fE, r ^F^(k) (56) 
Tin^^ TT JO k"^ 

a^{Za) 3 dx du dt Jt/2 , ft/2. 

Ttn-^ 2tt Jo I — x Jo Jo vl — ut aio^ a/o^ 

where 

/^^/2 = (1 - u)t[-3 + 5(1 - ut)] + 2x{l - ut) (57) 

and 

/3/2 = (1 - + ^(2 -t){l- Ut)] + X(l - t)[l + U{1 - t)]. 

Numerically we obtain 

2^E'hfs = 4.566(2) ^ ^' Ep. (58) 



4.2.2 Contribution to the Lamb Shift 

Consider now contribution induced by the spanning photon diagram with 
"left" self-energy insertion in Fig.2e to the Lamb shift. We will below repeat 
with minor changes considerations of the previous subsection. 

General expression for the spanning photon diagram coincides with the 
one in eg. (^9|) , only explicit expressions for the numerator structures slightly 
change. 

Unlike transformations in the previous sections we have to preserve tem- 
porarily small nonvanishing virtuality p = 1 — of external electron lines 
while combining denominators. This virtuality will be put to be equal to 
zero in the final formulae but it is necessary to preserve it on intermediate 
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stages to qualify spurious infrared divergences which appear in the subtrac- 
tion term below and cancel one another. Hence, we use instead of eg . (12^) 
slightly modified formulae 

A; = t[k\(l -nt) +t + ^ - + p]=tai (59) 

Explicit expression for the term Mp differs from the one used in the 
previous subsection only due to the change of free Lorentz indices 

Mp = l"{p + q)lo{p + q-k + l)7o(]3 + q+ l)7<x, (60) 

while we use another separation of different entries in the term A/"! here 

Afl = 2q{p + g)7o(p + q-k + l)7o(]3 + q + l)q (61) 

= (g^ + 2pg){[2g%o(p + q-k + 1)70 + 2gp7og7o] - 2gp7o^7o} 

+q'^2qfrfo{p + 1)70 + 4(pg)g]37o(]3 + 1)70 

= {q' + 2pq) [q'Afl, + Ml,] + q'Ml, + Ml,. 

Repeating calculations performed in the previous section we obtain (com- 
pare eq.(|5^)) 

A .3.^ _ . 2 Ml{q-Q)+Ml,{q-Q) 
txMUq~Q) , 'it{l-t)xMUq-Q) 



{q'-^if {q^-AiY 

m{l-u)eMUq-Q) 

{q^-^if ^' 

Performing next integration over momentum we obtain (compare eq.(|55D) 
for the electron factor 

r = 1 /' /' /' ,,,t(l-«)(3f-5) + x(t-3) 

16 JO 1 — X Jo Jo ai 

+—u[x{l - t) + (1 - u)t^{?,ut - 3m - 2) + — t(l -u){t^ + t- 6) 

ai 
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Next it is necessary to perform subtraction of the part proportional to 
from this expression. Note that the separate terms in the last brackets in 
eq. (|63|) would lead to infrared divergent integrals at k^ = if one omits small 
virtuality of the external electron line introduced in the beginning of this 
subsection. However, all integrals are perfectly convergent with nonvanising 
virtuality and we may use auxiliary identity from eq.( [L30| ) (for n = 2 and t] = 
1) to simplify the integral representation. After this simplification integral 
representation admits vanishing virtuality even at k^ = and we obtain after 
subtraction 



k4 k2 

1 dx /-i /-i u{l - ut) 



(64) 



' du C dti"^ — ^[x(3 - t) + (1 - u)t{h - 3t)] 
Jo aiam 



16 Jo 1 — X Jo Jo CLldlO 

1 I u 

2 ^ 2~ ' ~^ 5 
tti aio CLlCtiQ CLi 



-t^{l+t)u{l-u){l-ut){^ — + — ^) + -^[x{l-t) + {l-u)t^{-2-3u + 3ut)] 



:.u\l-u)t(l-t)}. 

af 

Substituting subtracted electron factor in eq.(||) and integrating over k 
we obtain 

2AEU.= ™ C^f^i^) fj^f^ (65) 

\m J nn-^ zvr Jo 1 — x Jo y/l — u Jo v 1 — ut 

'^ZO ^10 

where 

gl^^ = xu[l-t + 2{l-ut){3-t)]+ut{l-u)[t{-2-3u + 3ut)+2{l-ut){5-3t)], 

(66) 

gl^^ = -3u{l - u)t[2u{l -t)+ t{l - ut){l + t)]. 
Numerically we have 

2AEL„. = 5.056(1) (67) 
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5 Calculation of Contributions Induced by the 
Diagrams with Insertions of the "Right" 
Self-Energy Operator 

5.1 Diagrams Containing the Two-Loop Vertex 

5.1.1 General Expression for the "Left" Vertex with Right One- 
Loop Self-Energy Insertion 

Let us consider vertex diagram which is contained as a subgraph in Fig.2f. 
Insertion of the one-loop electron self-energy in the FY gauge corresponds to 
replacement (compare eq. (|23|) ) 

1 ? + g - ^) dx-^ — , , , , V, (68) 



p + q-k-1 47r Jo 1 - x + 2q{p - k) - 

Formal expression for the vertex has then the form 

d^q + N[/q 



K = 3(7^)' / dx 



1 X 



'477 Jo 1 — X J m"^ q^\q^ + 2gp](g^ + 2g(p — fc) + A;^ — 

where 



1-x' 



(69) 

Ni = riP + q + 1)7m(p + Q - kha, (70) 



N[ = 2q{p + q+ l)7^(p + q- k)q. (71) 
It is convenient to put down second numerator in the form 

N[ = {q' + 2pq)[q'N[, + N[,], (72) 

where 

nL = 27/., (73) 

N[, = 2^^{p-k)q. 

Combining denominators in eq. (|69D (compare eq.(p9|)) 

{l-t)q^ + t[{l-u){q^ + 2qp) (74) 

+u{q' + 2q{p -k)- -^)] = (g + Qf - A„ 

1 — X 
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Q — pt — kut, 

Ar = t[k^u{l -Ut)+t + ] = tUr, 

-L 2C 

we obtain representation for the bare vertex operator in the form 

Af = 3(f /' d.-^ r du /' dt I %2t^^^^^ (75) 

Nl,{q~Q) 2(1 -t)NUq-Q) 

where 

A^i = Ar{u = 1) = i[k^(l -t)+t + + p\= tari. (76) 

Next we perform shift of integration variable over momentum and eu- 
chdean rotation (and we shghtly change notation below as q further means 
euclidean shifted momentum) and obtain 

Af = C dx-^ C du [' dt r dq^q^i , , (77) 

" ^Att' Jo 1-xJo Jo Jo ^^^(g2 + A^)3 ^ ' 

[[3q^ - 2k''ut{l - ut) - 2(1 - - 2t{l - t)-t^{p - k - 1) 

+4(1 -t){l- ut)p^ip - ^ - 1) + 2(1 - t)\k] 

+ ^ ^ _ ^tjl-m^-l) 

^{q' + Ar,y^' {q' + Ar^y 

Momentum integration leads to 

A2 + A,i A2 ^ 2(l-t)(l-k2) ^ 

+2 log — 7^ + ^ ^7j. 

A^i + A^i a^i 
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Subtraction term is equal to 



Ai(0) =3(^)2 f'dx-^ f'du f'dt{2ti\og- 

^ An JQ 1 — X JO Jo 



ArO 2A2 + A,.o 

(79) 



+ 2[log — ^ , A ] + — }7m. 



QrO Ario A^ + A^IO OrlO 

where 

ttrO = O'r\]i=0^ 
Qrio = Ctrl|k=0- 

We extract explicit dependence on in eq.(^) after subtraction with 
the help of integration by parts 

, , ciro ,9 , t(l — u), 1 — ut. . . 

rfwlog— = kM du^ '-{u ), (80) 

ar Jo ar 

^ 1 ^'■10 -,2 dt t 
at log ■ 



/o ttri Jo a^i X + t{l — x) 

Then renormalized vertex with the right self-energy insertion acquires the 
form (compare eq. (|37|) ) 



H = 3(^)' /' dx-^ /' du /' dt{k'Ar7, (81) 
^ 47r JO 1 — X Jo Jo 

+i3r7^(p -k-l)+ CrP^,{p -k-l)+ Erlji}-, 



where 



2ut{l-t) , 2(l-nt)(^-t^) 2 

•^r \ 7 —f. TT, [P^j 

2t(l-t) 

Or , 

_ A{l-t)il-ut) 
c _ 2(l-t)^ _ 

Note that as well as in the case of the left self-energy insertion in eq.(P^) 
and eq.(^) all four functions Ar(k), Br{k), Cr(k) and Sr{k) are finite at 
k = 0. 
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5.1.2 Contribution to HFS Interval 

It is not difficult now to obtain contribution to HFS 

2AE{,^s = ^^^Ef{^) f'dx-^ I'du I'dt rrf|k|[A.(k) (83) 
vrn'^ TT^ Jo 1 — X Jo Jo Jo 

g,(k)-g,(0) 

where 

^^(k) - ^r{k)-£r{0) ^ ,r ^^(1-^^) ^ (1 - ut)[2u{l - t) - (1 - u)t'] 



(84) 

I }. 

ari{x + t — xt) 

Subtraction of £riS^) in eq.(p3D corresponds to subtraction of the contribution 
to HFS induced by the second order anomalous magnetic moment which is 
of lower order in parameter Za. 

After momentum integration we obtain three-dimensional integral for 
the contribution to HFS induced by both diagrams with nonsymmetric self- 
energy insertions in Fig.2f 

, a^iZa) , 3 , X du dt , . 

2AE{,ps = -^-1Ef{-3 + - / dx- / / -^== 85 

7m3 7T Jo I- X Jo ^/u Jo VI - Ut ^ ' 

}, 



ff ff 

r-ll/2 J3/2 



where 

///2 = -wt(l-t), (86) 

fl^ = (l-ut)[2u{l-t)-{l-u)t% 
Numerically we obtain 

2AE{,ps = -3A01{1)^^^Ef. (87) 
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5.1.3 Contribution to the Lamb Shift 

Contribution to the Lamb shift induced by the diagram in Fig.2f has the 
same form as in eg. p3|) , where 

Cf{k)-Cf{0) u{l-t){l-ut) 



(88) 



ar^rO arl[x + t{l — X)]'^ 

After momentum integration we obtain respective contribution to the 
Lamb shift in the form of the three-dimensional integral 

2AEl = m (!I^)3^!i^{-2 + dx^ (" du (" dt (89) 
m Txri'^ Ti Jo 1 — X Jo Jo 

where 

gi/, = -l+t, (90) 

g[/, = 2il-ut)iu-t^). 
Numerically we obtain 

2AE{ = -L017 1 m — f J ■ (91 



5.2 Spanning Photon Diagram with Symmetrical Self- 
Energy Insertion 

5.2.1 Contribution to HFS Interval 

Consider now spanning photon diagram with symmetrical self-energy inser- 
tion in Fig.2g. Unlike discussion of the vertex in the previous subsection 
we will consider here only the part of the spanning photon diagram rele- 
vant for hyperfine splitting. Insertion of the one-loop electron self-energy in 
the FY gauge corresponds to replacement described in eq. (^) in the central 
propagator in the spanning photon diagram. 

Formal expression for the spanning photon diagram has the form 
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/•A Nf, + Nl/q^ 

J m"^ + 2qp — p)^[g^ + 2g(p — k) + k"^ — — p] ' 

where 

A^F = 7"(]5 + g + 1)7m(p + ^ - ^)7^(P + g + 1)7-, (93) 
Nl = 2q{p + q + l)7^(p + q- k)^^{p + q + l)q (94) 
= 2(g2 + 2pg)^7^(g-p)7^. 

Note that we temporarily preserve in eq. (|92D nonvanishing virtuahty p = 
l—p^ of the external fermion lines. This virtuality will be put to be equal to 
zero in the final formulae but it is necessary to preserve it on intermediate 
stages to qualify spurious infrared divergences which cancel one another. The 
problem of infrared divergences is even more acute here than in discussion 
of the Lamb shift contribution induced by the diagram in Fig. 2e, since 
integral in eq.(^) contains more powers of integration momentum q in the 
denominator. We will see below that as a result separate parametric integrals 
giving contributions to the electron factor diverge (for vanishing virtuality) 
even in the case of nonvanishing exchanged momentum k and only total 
contribution to the electron factor turns out to be finite. 

Envisaging problems with infrared divergences it is convenient prior to 

integration to separate in the numerators A^^^. and Nl infrared safe terms Np^ 
and Nl^ 

m = N'pi + riP + 1)7m(p + q - khuiP + l)7<x (95) 
= Nl, + 47^(g - ^)7, = Nl, + Nl„ 

and 

Nl = q'Nl + 8{pq)\{q - k)^, (96) 
^q'Nl, + Nl,, 

Next we combine denominators as in eq. ([7^) but preserving nonvanishing 
virtuality of external electron lines so that explicit expression for A,, slightly 
changes and below 
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Ar = t[k'^u(l-ut)+t + + p]=tar. (97) 

1 — X 

After shift of integration variable q ^ q — Q (see eg. ([7^ ) we obtain 
retaining only those numerator structures which are relevant for contribution 
to HFS 



A(^)2 dx^— du2{l-u) f'dt3f (98) 

16 TT JO 1 — X JO Jo 

d'q ,NUq -Q) + [NU^l -Q) + Nliq - Q)] 



r(£q 



(g2 _ A,)4 
+4(l-t) ^^^^^-^H 



After integration over momentum we obtain 



+!fcM + H^(l^(3-< + kV()}. 



infrared unsafe terms in the central brackets in eq.(|99D may be easily 
exorcised with the help of the identity ineq. (|13CI|) . Substituting then expres- 



sion electron factor in eq.(|99D in eq.(|l]) and integrating over k one obtains 
contribution to HFS interval 



im-^ 2tt Jo 1 — X Jo Ju Jo J\ — ut 



a 



/2 



a/I — 

(100) 

where 

= 3 - 5Mt (101) 

and 
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/3% = -w(l-t)-t(l-Wt). 

Numerically we obtain 



AEf,ps = 2.682{1)^^^^Ef. (102) 
5.2.2 Contribution to the Lamb Shift 

Consider now contribution induced by the spanning photon diagram with 
symmetrical self-energy insertion in Fig.2g to the Lamb shift. We will closely 
follow considerations performed in the previous subsection. 

Formal expression for the spanning photon diagram has the same form 
as in eq.(^), the only difference is connected with the free Lorentz indices. 

We separate infrared safe terms in the numerators in the same way as in 
the previous subsection (see eq.(p5|) and eg . (|96|) ) and obtain repeating all by 
now standard steps 

\i^u{2eu - t - 1) + 2(3t - 2) 4(1 - t)(t - k^M^) 

~^ 2 ~^ 3 J ' 

This expression may be further simplified with the help of identity similar 
to the one in eq. (|130| ) but obtained with the help of integration by parts of 
the fraction t^(l —t)/ar instead of the fraction in eg. ( |12^ ). We then obtain 

a = - dx—^ C du{l - u) t dtt{ ^^'^^~^^ (104) 
8 Jo 1 — X Jo Jo 2ar 

, li''u{2t^U -t-1) 4^2(1 - t)2k2^ 



This expression is very convenient for subtraction of the low freguency 
asymptote since almost all terms in it contain factor explicitly. Performing 
subtraction we obtain 

U /:9(k) -£9(0) 3 /•! X /■! ,^ , /•! 



k4 k2 



/ dx—^ f du{l -u) I dtt (105) 
Jo 1 — X Jo Jo 



26 



3m(1 -Mt)(2t- 1) u{2t^u-t-l) Au\l-tf- 
2aaro al al 



Substituting now subtracted electron factor in eq.(|D and performing in- 
tegration over k we obtain contribution to the Lamb shift 

^ELmb=^[-^) (-) / dx—— / duMl-u) / dt 



m J nn'^ it Jo I — x Jo Jo a/1 — ut 

(106) 



where 



= 2-7t-3ut + 8M^^ (107) 

Numerically we obtain 

^Elmb = -0.14601(4) m (-^ j (108) 

6 Calculation of Contribution Induced by In- 
sertion of the Rainbow Self-Energy Opera- 
tor 

6.1 General Expression for the Rainbow Self-Energy 
Diagram 

Let us construct renormalized two-loop rainbow contribution to the electron 
self-energy operator in the FY gauge. To this end we take general expression 
for the one-loop bare self-energy operator in the FY gaugeQ 

S(P) = f / ^(^"^ + 2^)7..—^ J, (109) 



^We use dimensional momenta in this section 
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and perform substitution similar to the one in eg. (|68D , which describes 
insertion of the renormahzed self-energy operator in the FY gauge on the 
internal line. We then obtain 



Ei(p) = 3(iir|dx- 



X 



X 



it:- 



where 



(110) 



(111) 



2q{p + q)q = q\2q - p) + {2qpq + q^) = q^N^, + N^^. 



Combining denominators with the help of the Feynman parameters and 
performing shift of integration momentum^, Wick rotation and momentum 
integration we obtain explicit expression for the unrenormalized rainbow self- 
energy operator 



'3pH{p) 



3ppH'^{l-t) 
A 



where 



Hip) = log 



(112) 



(113) 



A, = t[-p'{l-t) + -^]. 

Term Ah is preserved on the background of the UV cutoff in eq. (|113|) due 
to the same reasons which were explained after eq.(^). It is convenient to 
rewrite H{p) in the form 

H{p) = H{m) + [H{p) - H{m)], (114) 

where first term is the momentum independent constant which in any 
case disappears after subtraction and the second term already admits limit 
of the infinite cutoff. We integrate the this second term over t by parts to 
get rid of logarithm and obtain 



^Note that Np + N^^ = ~3p does not contain linear in the integration momentum 
11, while integral depending on N^^ is 
under shift of the integration momentum. 



term, while integral depending on N^^ is UV finite, and, hence, no surface terms appear 
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sib) 



p dx 



X 



1 — X Jo 



dt 



Him) + 



pH\2-t) 
Ah 



(115) 



After two standard subtractions on the mass-shell followed by change of 
variables v = 1 — = v{l — x) and integration by parts over new variable v 
we come to the representation for the FY gauge renormalized rainbow mass- 
operator which is the most convenient one for subsequent calculation of the 
contributions to the energy splitting 



d^ p(l + + 2m 
-0^ m?-p'^i 



■ (116) 



Additional integration by parts over v leads to the representation in the 
form of one-dimensional integral 



a 



S"(p) = 9(-)^(p-m)M dv 
47r JO 



1 — V + log V 1 



p{l + v) + 2m 
m^ — p^v 



(117) 



= {p-mf{pap + ma^). 

Even remaining final integration may be performed analytically and one 
can obtain representation of the rainbow contribution to the electron self- 
energy in the closed form 

E«(p) = I^S^^(p) + ^(-)^^^^{2(p + m)r,{p) (118) 
Stt 16 vr 

-{p + 2m)J^2(p) - m^^^}, 

1 -p 

where p is the virtuality of the electron line (momentum p is arbitrary in 
this subsection) 



P 



1 1 
m — p 

m^ 



(119) 



S (p) is the one- loop electron self-energy (compare eq.( 



S^^(p) = -^(p-m)2r ^ 



P 



Air 



1-p (1-p)- 



■logp] 
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and functions Ti are defined as follows^ 

-[- 

6 



•^I(P) = -^\-- Li2(l - P) + p(logp - 1)], (120) 



•^2(p) = ^[y-Li2(l-p)]. 
6.1.1 Contribution to HFS Interval 

For final calculation of the contribution to the HFS interval we use represen- 
tation ([Tl6|) 



(121) 

c^\Za) ^ 9 f^dv 2 rd^ 1 + e _ 33 a\Za) ^ 



6.1.2 Contribution to the Lamb Shift 

With the help of representation (|116|) we obtain contribution to the Lamb 
shift interval 

AEl = m (^) ^ 3^ / T^{K(A;) + a^(A;)]-K(0) + a™(0)]} (122) 



m J nn^ \ vr^/Jo v Jo (1 — ^)^ 1 — + k^^ 

153 fmr\^a'^(Za)^ 
m ' ' 



40 \m J vrn^ 



^We use standard Euler dilogarithm Li2(p) here unlike function Li(p) in eq.(|lj) which 
was defined as the real part of the Euler dilogarithm. 
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7 Discussion of Results 



We presented above results of calculation of contributions to the Lamb shift 
and HFS of order cP'iZaf induced by all two-loop insertions in the electron 
line containing as a block one-loop electron self-energy graph. Calculations 
were performed in the FY gauge which is the most suitable one due to its in- 
frared smoothness. The formulae for different contributions obtained above 
admit numerical calculation with arbitrary accuracy. We consider it a bit 

premature to present here the sum of all contributions obtained above since 
considered set of graphs is not gauge invariant and comparison with the ex- 
perimental data have to await until all other contributions would be obtained. 
Nevertheless wc have chosen to present above calculations with sufficient de- 
tails both because of their volume and to present main technical tricks used 
in our work. 

Results of calculation of other contributions will be published separately. 

This paper was completed during the visit of one of the authors (M.E.) to 
the Penn State University. M.E. is deeply grateful to his colleagues at Penn 
State and especially to Prof. H. Grotch for kind hospitality. 

This work was supported by the Russian Foundation for Fundamental 
Research under grant #93-02-3853. Work by S.G. Karshenboim was also 
supported by a Soros Foundation grant awarded by the American Physical 
Society. 



Calculations performed in the main body of this paper were greatly facilitated 
by the infrared smoothness of the FY gauge. We would like to display in this 
appendix how cancellation of infrared divergences in the FY gauge occurs 
and derive a useful identity widely used above. 

Typical integral for the contribution to the electron factor with the worst 
infrared behavior has the form 



A 




(123) 
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where p is the on-shell electron momentum p = m and matrix A and 
function of the Feynman parameters u do not depend on q. 

Taking into account mass-shell condition for the vector p one easily ob- 
tains 



I = aA f dx... f dzf{x,...,z) I 

Jo Jo J TT^Z 



d'^q p 



2 



ttH q'^[q'^ + 2q{pri — kr) — uj]^~^^ 

(124) 

+ 2q{pr] - kr) - 

We omit below unessential for further considerations integration over the 
Feynman parameters and obtain after combining denominators, shift of in- 
tegration variable q ^ q — [pf] — kT)t and the Wick rotation 

/ = („+i)(-irv/;Ar/^{^^ (125) 



(g2 + A)"+2 

(g2 + A) ~ 



where 

A = {p^rft + k'^rH + u)t = ta. (126) 
Integrating over q one obtains 



1)" 2 , r 1 + (^/2)(1 - t) J(l - t)(mV + 't'^'^ 



'-^^ '-}. (127) 

Integration over t is in the general case infrared unsafe and may lead 
to infrared divergences if, e.g. variable u vanishes for some reason. Hence, 
it would be helpful to gain additional power of t in the numerator of the 
integrand. To this end it is useful to observe validity of the simple identity 

dt{l-t)_l-2t J{l-t)da _ l-2t ^t(l - t)(mV + '^^^^) 



dt a" a" 0*"+^ dt a"- a'^+i 

(128) 

Integrating the last term in eq. (|127|) with the help of the identity in 
eq. (|128|) one obtains 
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-rn' t^ll-Ulm- (129) 



n Jo 2 2 

We see that terms with minimal power of t in the numerator disappear 
when ^ = 2, i.e. in the FY gauge, thus making integration over t more 
smooth for low t. We used this trick with integration by parts abundantly in 
this paper and we put it down here for references in the FY gauge (i.e. when 

e = 2) 

/ dt{ 2n^ ' -} = S dt—. 130 

Jo a" a"+i ^ Jo a'^ 

Note that identity in eg. ( P-30|) is only a representative of a large family 
of identities which may be obtained in the same way but choosing different 
numerators in eq. 
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Figure Captions 

Fig.l. Six gauge-invariant sets of diagrams producing contributions of order 
a\Zaf to HFS and Lamb shift. 

Fig. 2. All graphs with two radiative photons producing contributions of order 
o?{Zaf to HFS and Lamb shift. 
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